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, Abstract: QCD perturbation theory using quark and gluon propagators that have an 

\ extra term at zero four-momentum is equivalent to standard PQCD in the presence of 

\ massless quark and gluon pairs in the perturbative vacuum. We verify at low orders in 

' perturbation theory that {il^ip) 7^ and that quarks acquire a constituent mass. Thus chiral 

Q_i' symmetry is spontaneously broken, and the conservation of the axial vector current ensures 

a massless pion. The gluon pairs generate an {F^^F^^'^) expectation value and a tachyonic 
gluon mass. Lorentz and gauge symmetry 
PQCD is unaffected. The modified perturl 
of QCD at both short and long distances. 



D . gluon mass. Lorentz and gauge symmetry is preserved, and the short-distance structure of 

PQCD is unaffected. The modified perturbative expansion thus contains relevant elements 
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1. Perturbing in the presence of parton pairs 



We wish to call attention to the interesting properties of a Perturbative QCD (PQCD) 
expansion where the free, massless quark and gluon propagators have an extra term at 
vanishing four-momentum. In Feynman gauge, 



+C,{2.f5\p) 



' +C,{27,)H\p) 



p^ + ie 



(1.1) 



[1.2) 



A corresponding modification of the ghost propagator is irrelevant since the ghost-gluon 
vertex is proportional to the ghost momentum. 

The gluon propagator modification (|1.2[) was proposed in [l|, motivated by a Lorentz 
and gauge invariant formulation of QCD vacuum effects. In later work it was shown that 
such a modification would arise if gluon pairs are present in the asymptotic wave function 
Q , and some of its phenomenological consequences were studied Q . Independently, it was 
noted in Q that perturbative expansions using different ie prescriptions at the poles of the 
free gluon propagator are equivalent to standard expansions with gluons added to the in- 
and out-states. This was inspired by earlier studies of the relation between modifications 
of the on-shell propagators and boundary conditions. 

It should be emphasized that a perturbative expansion using propagators such as ( |1 . 1| ) 
and ( |1.2| ) is formally as justified as the standard expansion. Standard perturbation theory 
uses empty, perturbative vacua as in- and out-states at asymptotic times {t = iLcxd). In 
the infinite time separating the empty states from the scattering event they will relax to 
the true QCD ground state (assuming the overlap to be non- vanishing) . At finite orders 
of PQCD one sees, however, only the beginning of this relaxation. Thus if the true ground 
state differs essentially from the perturbative vacuum, e.g., by containing quark and gluon 
condensates |^, the perturbative expansion will not capture this physics. This is a likely 
reason for the failure of standard PQCD at long distances, where soft condensate partons 
have important effects. 

It is thus natural to consider PQCD with asymptotic states that contain (free) quark 
and gluon pairs. Assuming again an overlap of those states with the true ground state the 
resulting expansion will be a priori as justified as standard PQCD. The addition of quark 
and gluon pairs to the asymptotic states implies no change in the QCD lagrangian nor in 
the short distance properties of the theory (including renormalizability) . Such asymptotic 
states may on the other hand give rise to non- vanishing quark and gluon condensates, 
'constituent' parton masses and a spontaneous breaking of chiral symmetry. 

There are phenomenological indications that perturbative methods are relevant for 
long distance QCD physics. Constituent quarks are effective degrees of freedom in hadrons. 
The transition from short distance (quark and gluon) to long-distance (hadron) physics is 
smooth 0, 0. This suggests that the strong coupling freezes, UsiQ"^ = 0) ~ 0.5, and that 
confinement and chiral symmetry breaking is due to quark and gluon 'condensates' in the 
QCD vacuum [0, rather than to a strong coupling regime. 
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The abstract contains a detailed study of how gaussian wave functions of fermion 
and scalar fields at asymptotic times t = ±T — > ±00 imply free propagators of the form 
( |1.1| ) and ([L.2|). The coefficient of the zero-momentum contribution is determined by the 
gaussian, and vanishes when the wave function is that of the ground (empty) state. In 
general the gaussian may be expanded to reveal an indeterminate number of particle pairs. 



2. Vacuum expectation values and constituent masses 

The consequences of the modifications ( |1.1| ) and ( |1.2| ) to the quark and gluon propaga- 
tors are profound. Due to the infinite number of pairs the modified expansion is non- 
perturbatively related to standard PQCD. 

We find that the propagator modifications give rise to non-vanishing {ipip) and {F^yF^'^) 
expectation values. Thus chiral symmetry is spontaneously broken. Quarks gain a con- 
stituent mass whereas gluons become tachyonic due to interactions with the condensates. 
We verify (at lowest non-trivial order) that the altered boundary conditions do not affect 
the conservation of the vector and axial currents, d^J^ = d^J^ = 0, which follows from 
the symmetries of the lagrangian. Thus the pion is massless, despite the masses of its 
constituent quarks. 



2.1 Quark and gluon condensates 



The quark propagator ( |1.1| ) leads to spontaneous chiral symmetry breaking already at 
zeroth order in perturbation theory. 



(V(x)V'(x))=5^ 



5 



a, A 



exp 



(2^ 



r(p)5rc^(p) 



(2.1) 



where is the number of colors. In effect, the divergently large number oc 1/m of zero- 
momentum quark pairs implied by Eq. ( A.23| ) compensates the factor m in the numerator 
of the propagator, giving a finite result in the m — > limit. Thus chiral symmetry is 
spontaneously broken by the asymptotic state. 

In the case of the {F^yF^'^) expectation value the O ((7") contribution vanishes since 
the free part of F^y = d^Ay — dyA^ — gfahc-^^^-^u is proportional to the gluon momentum. 
Furthermore, at O ((7^) the gluon loop integrals vanish (when dimensionally regularized) 
due to the absence of a mass or momentum scale. This leaves only the two-loop diagram 
involving a 4-gluon coupling (Fig. la), with the Cg part of the propagator taken in both 
loops. The result is pi, 0] 



{F,y{x)F^^%x)) 



5^ 



12g'^N{N'^ - 1)C^ 



j{-p)Dgj{p) 



(2.2) 



The gluon condensate thus arises specifically from the non-abelian 4-gluon interaction. 
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Figure 1: (a) The only diagram which contributes to (F^i^F^^'^) at O (g^)- The gluon condensate 
term cx Cg in ( |l.2| ) is taken in botli loops, (b) Quark propagator correction. (c,d) Gluon propagator 
corrections. 



2.2 Quark constituent mass 

The quark propagator correction S(p) shown in Fig. lb gets a contribution from both the 
quark condensate term Cq in (1.1) {k = 0) and from the gluon propagator modification 
Cg of ( [1.2| ) (k = p). Only the former generates a proper mass term, making TiSq{p) ^ 0, 
where (for p ^ 0) 

It is straightforward to write down the quark condensate contribution to the diagram of 
Fig. lb, 

(p) = ^9'C^Cq4f- (2.4) 

^ p'^ + l£ 

where Cp = {N'^ — 1)/2N. As expected, it is suppressed in the short distance (p^ —f oo) 
limit. The pole in the propagator (2.3) is at the 'constituent' quark mass p^ = M^, with 

Mq = (Vc^lClj (2.5) 
The gluon condensate correction to the quark propagator is similarly |p, 

Ef^ip) = -2g^CpCg^^6^^ (2.6) 
" p'' + le 
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and influences the quark constituent mass Mq, which is now is obtained as the solution of 
the cubic equation, 

- Ig^CpCgMg = ±Ag^CFCg (2.7) 

2.3 A tachyonic gluon 

The gluon propagator corrections shown in Figs. lc,d were evaluated in Ref. [^. The 
sum of the two diagrams has the transverse structure required by gauge invariance. The 
resummed propagator thus has the form (for p ^ 0) 

Df'^^ip) = ,s^^Z£Z±p!^ _ ,5"^^ (2.8) 

where the gluon condensate contribution is 

Iig{p) = 2g^NCg (2.9) 

There is no contribution to the gluon propagator from the quark condensate term Cq at 
O (5^) . The insertion of this term in the quark loop correction gives a vanishing trace due 
to an odd number of 7 matrices. 

A modification of the ghost propagator at p = also does not contribute to II^ since 
the ghost-gluon vertex vanishes at zero ghost momentum. In Ref. where asymptotic 
states of gluons with p 7^ were considered, the ghost contribution was non- vanishing and 
essential for gauge invariance. 

As was found also in |8|, the gluon propagator has a pole at negative p^, corresponding 
to a 'tachyonic' gluon mass, 

M2 = -2g^NCg (2.10) 

provided Cp > 0. This is likely to have important consequences. A quark or gluon 
can acquire an absorptive part by emitting a spacelike gluon with p^ ~ < 0. Thus 
colored partons cannot propagate long distances due to their interactions with the quarks 
and gluons in the asymptotic state. Any stable eigenmodes of propagation are likely to be 
color singlets, which do not emit soft gluons. 





(a) (b) 

Figure 2: Corrections to the color singlet (axial) vector current. A third diagram where the gluon 
couples to the antiquark is not shown. 
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2.4 Conservation of vector and axial vector currents 

The breaking of chiral symmetry ( p.lD resulted from our modification of the asymptotic 
states. The QCD lagrangian is unchanged and remains chirally symmetric. Consequently 
the symmetry is broken spontaneously and there must be a massless goldstone boson, the 
pion, which is annihilated by the axial current, 

{0\4{x)\n{p)) = -ip^'Uexpi-zp-x) (2.11) 

The symmetry of the lagrangian implies d^J^ = 0, i.e., m-,^ = 0. Current conservation 
should hold order by order in the modified perturbative expansion. We shall verify it for 
the diagrams of Fig. 2. 

The quark condensate term Cq of ( [L.l| ) must be used in two of the four quark prop- 
agators in Fig. 2 to have an even number of 7 matrices in the quark loop. Consequently 
the two-loop integral reduces to a number, and we find, for a vector {—iej^) coupling of 
the current, 

(p) = 16e^g\N^ - 1)^ (5^^^ - ^) (2.12) 

The transverse structure ensures vector current conservation and gauge invariance. The 
result for an axial vector current coupling (—^6757'^) is the same, up to a sign. 

The contribution of the gluon propagator modification (|l.2| ) on the gluon line in Fig. 
2 reduces the two loop integrals to a single (regular) integral. For a vector coupling we find 

^r(^') = ^(^' - 1)^9 (5^' - ^) (2-13) 

We have again Pf^Ag'^{p) = as required by current conservation. The result is the same 
for an axial vector current. 

3. Concluding remarks 

The perturbative expansion is a powerful tool in field theory. It is specified by the la- 
grangian and by the state which one chooses to expand around. The sum is formally 
independent of this state provided it is taken at asymptotic times {t — > ±00) and has 
an overlap with the true vacuum. Hence it is interesting to consider the properties of 
perturbative QCD expansions with non-empty asymptotic states. 

A PQCD expansion can be phenomenologically useful provided quarks and gluons are 
relevant degrees of freedom for the physics being studied. In the short distance regime 
this is clearly so, and PQCD has been impressively successful. While hadrons in many 
respects are good degrees of freedom at long distances, there are indications that partons 
maintain their identity even in this regime. Thus hadrons can be characterized as loosely 
bound states of constituent quarks, which differ from current quarks only in their mass. It 
has also been noted that the transition from the short distance to the 'confining' regime 
is smooth 0. This is consistent with a picture where the strong coupling freezes at 
a perturbative value ag^Q"^ = 0) ~ 0.5 and the novel long distance effects result from 
interactions with quark and gluon condensates in the physical ground state. 
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We previously studied |^ |^ the effects of including finite momentum quarks and gluons 
in the asymptotic states. Such states are not boost invariant and Lorentz invariance is 
thus implicit, i.e., the convergence of the expansion is frame dependent. This may not be 
a disadvantage when a choice of frame is anyway implied, as in the calculation of bound 
states wave functions. However, the possibility of retaining explicit Lorentz invariance by 
including only massless partons with vanishing momenta [||, ^ is attractive. This is also 
similar to the framework of the QCD sum rules Q, where no momentum is exchanged with 
the vacuum condensates. 

The modification of the gluon propagator ( |l.2D in Feynman gauge implies a corre- 
sponding modification of the ghost propagator. This was shown to give a gauge invariant 
result in a one-loop calculation of the static potential Q. However, since the ghost gluon 
vertex is proportional to the ghost momentum a modification of the ghost propagator at 
p = is irrelevant. The results obtained here and in Ref. Q without ghost modifications 
are indeed consistent with gauge invariance. 

In this paper we extended the considerations of Refs. [Q, ^ to quarks and found 
that chiral symmetry breaks spontaneously. We made non-trivial checks that identities 
following from symmetries of the lagrangian are respected by the condensate contributions. 
In particular, we saw how the pion can stay massless despite the finite mass of its constituent 
quarks. The quark mass contribution (Fig. 2a) is exactly cancelled by the interaction 
between the quarks in (Fig. 2b). 

The fact that the gluon propagator acquires a pole for < (a tachyonic gluon mass) 
|3|, P is likely to have important consequences. It implies that quarks and gluons can 
propagate only over a limited range in the condensate, as the emission of long wavelength 
spacelike gluons generates an absorptive part. Since the colored condensate will be more 
transparent to color singlet bound states it is not farfetched to speculate that hadrons are 
the only stable eigenmodes of propagation. 
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APPENDIX 

A. Influence of the asymptotic state on the free propagator 

The generating functional Z[J] of Green functions in (scalar) field theory can be expressed 
in the form 

1"Jint 



Zint[J\ = exp 



Z[J] (A.l) 



where J is the source of the field Sinti^p) is the interaction part of the action and the 
free functional is explicitly given by the free propagator D{p), 

Z[J] = J V[(j)]exp!^i j d'^x[Coi(j)) + Jix)(l)ix)]^ (A.2) 
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exp 



(2^ 



J{-p)D{p)J{p) 



(A.3) 



The boundary condition at = itcx) of the path integral ( A. 2 ) is usually specified 
implicitly by tilting the integration path to the real axis, giving the Feynman propagator 



Df{p) 



p 



2 — m? + ie 



(A.4) 



In the operator formalism this propagator is obtained by requiring the perturbative vacuum 
to be empty, a|0) = 0. 

We wish to show how the inclusion of parton pairs in the asymptotic state modifies 
the prescription at the p^ = m? pole of the free propagator. We consider first the case of a 
scalar field and then that of a fermion, making use of the results in The generalization to 
the free gluon and quark propagators of QCD is straightforward. The full QCD generating 
functional is given by the free one as in Eq. (A.l). The standard Feynman rules of PQCD 
are thus modified only through the change in the free propagators. 

A.l Free scalar field 



The space part of the gaussian path integral in ( |A.2| ) is diagonalized in 3-momentum 
space. The following analysis is thus independent of p. For simplicity, and since we shall 
be interested in zero momentum pairs, we assume p = and treat only the time/energy 
dimension explicitly. 

To see the influence of the asymptotic states on the free propagator we need to first 
calculate the path integral over a finite time interval, —T < t <T, and then take T — > oo. 
We define to be 



Z[J] = / V[<p{t)] exp 



dt<j){t) 



+ m 
_ 

2 



(j){t)+i / dtJ{t)(t){t) 



T 



cp\-T)+cp\T) 



(A.5) 



Here Cg > parametrizes the gaussian wave functions at t = itT. Since the free scalar la- 
grangian describes a set of harmonic oscillators we expect that the ground state is obtained 
for Cs = m. We find below that precisely this value of Cg gives the Feynman propagator. 
Other values of Cg correspond to a superposition of wave functions containing particle 
pairs. 

Describing the field through its Fourier components. 



1 



oo 

2T ^ 

p=—oc 



, exp(— ivrpt/T) 



(A.6) 



implies a periodic boundary condition, 

m = ci>{-T) 



2T 



(A.7) 
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Introducing this relation via a (^-function constraint gives (up to an overall normalization), 



Z[J] 



da d4>{T) Yl d(t)p exp <j ^ J]] (p-p 



p=— oo 



-lfa]cl)p-Cscl)\T)+iacP{T) 



(Ai 



The result of the gaussian integration is 



Z [J] = exp ■ 



oo J J 

E'-' ~p^p 

4t— — 7T7 

p=-oo 



1 ^ 1 1 [^Y-J 



(A.9) 



This is a well-defined expression if mT ^ np for any integer p. In order for the T ^ oo 
limit to be smooth we take 

mT= I A^ + ijvr (A.IO) 



which ensures that the denominators vanish in the middle of a discrete momentum interval, 
and then let ^ oo. 
Using 

E°° I ( — 1)^ vrcosfTra) 

— = cos(7rr^ ^ — 
— 

p=— oo 

implies with ( A.1C| ) that 

y — = 



p"^ — a"^ a sin(7ra) 



(A.ll) 



(A.12) 



p= — oo J'2 

in Eq. ([A. 91) . In the continuum (N oo) limit the source Jp is slowly varying compared 
to (—1)''' in the last term of ( A.Dl) . Hence that term only contributes for p ~ ttiT/tt, where 
the denominator is rapidly varying. Thus, using again ( A.ll[) , 



lim — - — 



-l)PJn 



2p2 2 



j(^) (_1)P 



2T 



V _ ™2 
p '"' 



J{m) 
2m 



(-1) 



N+l 



(A.13) 



where according to (A. 6) we indicate the continuum energy E as Jp = J{E = 'Kp/T). 
Similarly we note that the continuum limit of the first term in ( A.9 ) is 



lim 

N^oo AT 



V — ^ 

Z-^ TT^P' 
p= — oo J'2 



J—pJp 



dE J{-E)J{E) 
2^ E"^ -m? 



(A.14) 



. J{—m)J{m) 



oo ^ 

V - 

Z—d ■k'^p'^ 



p=— oo 7^2 



where V denotes the principal value. The second term represents the contribution close to 
the pole of the propagator where the sum may not be approximated by the integral. This 
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term vanishes due to ( A.lCj ) and ( [A.ll[) . Collecting the above results gives 



Z[J] = exp<^ 



dE J{-E)J{E) _ _Cs 



Am? 



exp 



dE 
2^ 



J{-E) 



+ 



1 



_ _|_ ^£ 2m 



9l 
m 



1 2tt5{E 



m 



(A.15) 

J{E) 



Note that J[—m) = J(?ti), since only the combination Jp + J_p enters in the pole contri- 
bution ( [Al3D . This is a consequence of the periodic boundary condition ( [A.71 ) . In general 
is related to <j){—t) through (pp — > (p-p, hence the asymptotic momentum space field 
(p{m) is real. 

Eq. ( A.15| ) shows that the Feynman propagator is modified for Cg m, analogously 
to the gluon propagator in Eq. (|1.2|). As we noted above, Cg = m corresponds to imposing 
the ground state wave function as boundary condition at t = ±T ^ zboo. For Cs ^ m the 
gaussian wave function may be expanded in powers of (C^ — m)4?{-^T) and interpreted as 
a state with additional particle pairs. 

It is straightforward to extend this analysis to fields (j)x{t,p) having non- vanishing 3- 
momenta p and several spin components A. If the argument of the gaussian wave function 
is of the form (pxi^T, —p)(j)xiiT,p) it will represent spin zero pairs with vanishing total 
3-momentum. 

An application to the gluon propagator is more subtle due to gauge invariance. 
In the limit of vanishing mass m we must let Cg — > m in ( A.15| ). The sign of Cg in ( |1.2| ) 
then depends on whether Cg approaches m from below or above. In the QCD analysis of 
Ref. m it was found that Cg > 0. We use Feynman gauge as in 1^, Gauge independence 
was explicitly verified |^] for Feynman and Coulomb gauge in a one-loop calculation of the 
static potential with boundary conditions modified at p ^ 0. In that case a modification 
of the ghost propagator analogous to was essential in Feynman gauge. It is, however, 
irrelevant in our present (p = 0) analysis since the ghost-gluon vertex is proportional to 
the ghost momentum. 

A. 2 Spin i field 

We take the finite time, free generating functional analogous to ( |A.5| ) for a spin | field to 
be 



Z[C,C]= / V[i;,^P]exp 



dt 



- C 



^(_r)v(-r) + v(r)^(r) 



(A.16) 



where the fields ^ and sources C, C are Grassmann valued Dirac spinors, and the wave 
functions at t = itT are parametrized by the constant Cj. For notational simplicity we do 
not indicate the 3-momentum (and take p = 0) . The Dirac algebra is diagonalized using 
fields Xo-,A5 Xa\ where a = + (— ) is an (anti) particle label and A =t,i indicates spin 
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component, 



(A.17) 



Here n^(T) = (1, 0, 0, 0), u^{i) = (0, 1,0, 0), and similarly in terms of the lower components 
for v{X). 

Defining the source 770- x correspondingly. 



C(t)=7°Eh(^)^+-A(i)+t;(A)r?_,A(t)" 



(A.18) 



the free fermion functional ( A.16| ) becomes 



Z[r],r]] = / D[x,x]exp <^ ^ 



cr,A 



dt 



- C/E^[^-,A(-r)x<x,A(-T)+x.,A(r)x<x,A(T)]| (A.19) 
In terms of the Fourier components defined as in ( |A.6| ) we get, 

Z[ri,r]] = j Yl (iXp dXp exp - ^ Xp^pgXg + ^"^p^p + ^p'^p) 



det(^) exp 



P.9 



(A.20) 



where we suppressed (the sum over) a, A in the exponent. The inverse of the matrix 

aC 



A 



pq 



(A.21) 



2T \T 

can be readily guessed from the corresponding result ( [A .91 ) for scalar fields. Imposing again 
( A.ICP and using ( A.ll| ) we find (for each a, A and ignoring the overall normalization) 



Z[ri,r]] = exp 



- y 



VpVp 



+ 



Vpr]q 



p=—oo 



^y^(_i)p+9 

p — amT/iT ' 27r2 ^—^ {p — amT /iT){q — amT /tt) 



P,Q 



(A.22) 



In the continuum (N — > 00) limit we get, similarly to the scalar case above. 



Z[r],r]] = exp<^ ^ 



iV 



2iT E — am 



+ 



(A.23) 



exp 



00 JT71 _ 

— C{E){E^^ + m) 



V 



E^ 



+ CfTT5{E'^ - rr? 



CiE) 



where Cf = 1 gives the Feynman propagator. In the m ^ limit C(ibm) C(0) the 
£^7^ term cancels in the sum of the E = itm pole contributions. Thus 

Z[rj,7jr=%^plv r '^l{E)'-^Q{E) + \Cfmm] (A.24) 



which gives a free fermion propagator of the form ( |1 . Ij) . 
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